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A small value of the spin gap in quantum antiferromagnets with strong frustration makes them susceptible 
to nominally small deviations from the ideal Heisenberg model. One of such perturbations, the anisotropic 
Dzyaloshinskii-Moriya interaction, is an important perturbation for the S = 1/2 kagome antiferromagnet, one 
("■^ , of the current candidates for a quantum-disordered ground state. We study the influence of the DM term in a 

^vj ■ related one-dimensional system, the sawtooth chain that has valence-bond order in its ground state. Through 

a combination of analytical and numerical methods, we show that a relatively weak DM coupling, 0.11 5 J, is 
sufficient to destroy the valence-bond order, close the spin gap, and turn the system into a Luttinger liquid with 
algebraic spin correlations. A similar mechanism may be at work in the kagome antiferromagnet. 
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I. INTRODUCTION 



Antiferromagnets with S = 1/2 and on non-bipartite lattices are considered viable candidates for exotic ground states and 
1* , excitations. Geometrical frustration and strong quantum fluctuations tend to suppress long-range magnetic order. The resulting 
5— I ■ ground state does not break the symmetry of global spin rotations, but its exact properties remain subject of vigorous debate, with 
c/3 \ proposals ranging from valence-bond crystals that break some lattice symmetries'"-' to valence-bond liquids that fully preserve 
the symmetry of the Hamiltonian."'"^ A spin-liquid state with an energy gap to all excitations may further possess a hidden 
2 ' topological order. Several antiferromagnetic materials without long-range magnetic order well below the characteristic Curie- 
Weiss temperature scale have been discovered recently, most notably herbertsmithite Cu3Zn(OH)gCl2,** where no magnetic 



f, ' order has been detected down to 50 mK,^ '-^ even though the exchange interaction is estimated to be J = 180 K. The material 



is a "structurally perfect"'''''* realization of the S = 1/2 Heisenberg antiferromagnet on kagome, a network of corner-sharing 

Q ' triangles. Fig. 1(a). 

Q [ While most of the theoretical studies of quantum antiferromagnets deal with the pure Heisenberg model with nearest-neighbor 

'~~' • exchange, real systems inevitably deviate from this idealization. Frustrated magnets in particular are sensitive to various nomi- 

_^ [ nally weak perturbations. In this paper, we deal with the Dzyaloshinskii-Moryia (DM) interaction,'^'^ the antisymmetric version 

h>. ■ of the Heisenberg exchange induced by the spin-orbit coupling. The Hamiltonian of such a system is 

0^ . H = J2[J S. • S, + D,, • (S, X Sj)]. (1) 

?i : '' 

f — . In herbertsmithite, the DM term is allowed by the crystal symmetry. The in-plane and out-of-plane components of the DM vector 

l^ ' Dy on kagome are shown in Fig. 1(b) and (c). From eSR measurements,'^ the DM vector has the magnitude D = 0.08 J and is 

^-H , dominated by the out-of-plane component, whereas the in-plane component is small. Din = O.OIJ ± 0.02 J. The DM term can 

^~~^ ' be gauged away by an appropriate rotation of the local spin axes,'^'^ provided that its "line integral" vanishes for any closed 

loop abc . . . yza: 



Bat + 'Dbc + --- + T>y, + D,„ = 0. (2) 

It can be seen from Fig. 1(b) that the in-plane component satisfies Eq. (2) and thus can be gauged away. The out-of-plane 
component cannot be removed in this way and thus represents a physical perturbation. In this work, we concentrate on the 
out-of-plane component of D. 

A growing evidence from numerical studies^^"^'' indicates that the pure Heisenberg model, D — Q, has a S* = ground state 
with a small but finite energy gap for S* = 1 excitations, with estimates ranging from A = 0.05 J to 0.15 J. These values are 
comparable to the strength of the DM term, so it is plausible that the low-energy properties of herbertsmithite are influenced by 
the DM interaction. 

The effects of the DM interaction on the kagome antiferromagnet were first studied by Rigol and Singh^^-^^ in order to explain 
low-temperature paramagnetism in herbertsmithite: an upturn in magnetic susceptibility at low temperatures^' seems to indicate 
the absence of a spin gap. Tovar et alr^ concluded that a finite DM term could be responsible for the non-zero susceptibility 
observed in experiment even if the spin gap remains finite. A study employing exact diagonalization"^ showed that a sufficiently 
strong DM term, D > Dc~ 0.10 J, induces long-range magnetic order in the ground state, with magnetic moments lying in the 
plane. This was later confirmed by employing the Schwinger-boson approach. ^''•^' The ordering tendency is easy to understand 




FIG. 1: (a) Kagome lattice, (b) and (c) In-plane and out-of -plane components of the DM vector Dij shown for directed links (i -^ j) on 
kagome. 



by turning to the classical variant of the Heisenberg model. There, the out-of-plane D vectors shown in Fig. 1 (c) lift the extensive 
degeneracy of the classical ground states leaving a q = ground state that spontaneously breaks the remaining 0(2) symmetry 
of the DM Hamiltonian (1). Later numerical work^^ turned up some evidence that the system may have an intermediate phase 
between Dd ~ 0.05 J and Dc2 ~ O.IOJ, where 5*2 = 1 excitations become gapless but the spin 0(2) symmetry remains 
intact. In the absence of an obvious order parameter that would uniquely identify the intermediate phase, the authors of Ref. 32 
concluded that the appearance of an intermediate phase might be a finite-size effect. Further work in this direction is required to 
elucidate the nature — and even the existence — of the intermediate phase and its possible relevance to herbertsmithite. 

In our previous work,^^ we have shown that the 5* = 1/2 Heisenberg antiferromagnet on kagome can be viewed as a collection 
of fermionic spinons — topological defects with S = 1/2 — moving in an otherwise inert vacuum of valence bonds. The spinons 
interact with an emerging compact U(l) gauge field whose quantized electric flux is related to the valence-bond configuration 
through Elser's arrow representation.'''* Spinons carry one unit of the U(l) charge against a negatively charged background. 
These features are reminiscent of the picture of fermionic spinons proposed earlier by Marston et a/. ''^^ and Hastings'*, who 
used the Abrikosov-fermion representation for spin operators. It is worth pointing out that the Fermi statistics of spinons is not 
postulated ad hoc in our approach but rather arises naturally as the Berry phase of valence bonds that are moved in th e process 
of spinon exchange. We have further shown that strong, exchange-mediated attraction binds spinons into small and heavy 5 = 
pairs and that low-energy S* = 1 excitations result from breaking up a pair into "free" spinons. Thus the spin gap is determined 
mostly by the binding energy of a pair, which we estimated to be 0.06 J. 

From this perspective, one potential route to the closing of the spin gap could be via the destruction of the two-spinon bound 
state in the presence of a sufficiently strong DM term. That, however, appears unlikely for two reasons. First, the factors setting 
the pair binding energy — the spinon hopping amplitude and the strength of exchange-mediated attraction-are both of order J, 
so it is hard to see how a fairly weak coupling D ~ 0.05 J to O.IOJ can disrupt the pairing. Second, a quantum phase transition 
to a state with long-range magnetic order can be viewed as Bose condensation of magnons,^^ quasiparticles with Sz = I and 
there are no low-energy excitations of this kind in the pure Heisenberg model. Although one could think of condensing pairs 
of spinons with Sz = 1, this route runs into another difficulty: such an object would carry a double U(l) charge, whereas a 
magnon is expected to be neutral. Put simply, a pair of spinons is a topological defect whose motion affects the valence-bond 
background, which is uncharacteristic of magnon motion. 

A possible way out is to postulate that the condensing objects are pairs consisting of a spinon and its antiparticle. Such a 
composite object would have zero U(l) charge and be topologically trivial, like a magnon. In the pure Heisenberg model, the 




FIG. 2: (a) The sawtooth chain, (b) and (c) Its valence-bond ground states, (d) Spin- 1/2 excitations: kink (left) and antikink (right), (e) 
Orientation of the DM vectors Dij. (f) The ground state of the classical model has a commensurate magnetic order with the wavenumber 
<j/27r = -1/3. 



energy cost of creating a spinon and its antiparticle is approximately 0.25 J. ^^ As we will see, the DM term lowers the kinetic 
energy of both spinons and their antiparticles. It is thus reasonable to expect that, at some critical coupling strength Dc, the 
energy cost of adding a pair vanishes. 



To test this scenario, we have studied a toy version of the kagome antiferromagnet known as the sawtooth spin chain,^'^'^^ a 
one-dimensional lattice of corner-sharing triangles. Fig. 2(a). To make a connection with kagome, exchange couplings are set 
equal for all bonds. At D = 0, the chain has two valence-bond ground states. Fig. 2(b) and (c), that spontaneously break the 
mirror reflection symmetry. Spin excitations are topological defects: domain walls with spin S = 1/2, Fig. 2(d). The domain 
walls come in two flavors: kinks have zero energy and are localized, whereas antikinks are mobile and have a minimum energy of 
0.215 J. ""* These excitations can only be created in pairs by a local perturbation acting in the bulk. As we discussed elsewhere,"**^ 
spinons of the kagome antiferromagnet have similar properties, with one notable exception: the ground state of the sawtooth 
chain is free from the defects, whereas kagome has a finit e concentration of antikinks (1/3 per site) bound into S — pairs. 

We have studied the sawtooth spin chain with exchange and Dzyaloshinskii-Moriya interactions, Eq. (1). The Dij vectors 
had the same length and a uniform out-of-plane orientation preserving the translational symmetry of the chain as shown in 
Fig. 2(e). Qualitatively similar results were obtained for the staggered choice of D^j, but we will not provide the details here. 
The introduction of the DM term preserves the mirror symmetry of the Hamiltonian (it inverts the x coordinate of the lattice and 
the Sy and 5*2 components of the spins), so that the notion of a valence-bond order that spontaneously breaks this symmetry is 
still valid. The valence-bond order survives to a finite value of the DM coupling. 

As described below, kinks become mobile in the presence of a DM term. Their minimal energy becomes negative, growing 
linearly with D. The minimal energy of an antikink remains unchanged to the first order in D, so one can expect that the 
minimum energy of a kink-antikink pair will vanish when D reaches a critical value Dc of the order of the initial spin gap, 
0.215 J. In Sec. II, we describe a calculation of the spinon spectrum in the presence of a nonzero D, from which we obtained an 
estimate of the critical DM strength, Dc — 0.087 J. For D > Dc, spontaneous creation of kink-antikink pairs leads to a finite 
concentration of topological defects, which obliterates the valence-bond order and restores the reflection symmetry of the lattice. 
This scenario is reminiscent of quantum phase transition at the end of magnetization plateaus in the S* = 1/2 Ising-Heisenberg 
chain**' and in a frustrated two-leg ladder.'*^. In both of those models, the condensation of domain walls turns a state with a 
broken translational symmetry and gapped excitations into a gapless phase with incommensurate spin correlations decaying as a 
power of the distance. Exact diagonalization calculation for the sawtooth chain with DM interactions, described in Sec. Ill, are 
consistent with this scenario. 



II. SPINON DISPERSIONS 

A. D = 

We briefly review the physics of the sawtooth chain in the pure Heisenberg model without the DM term.''^^° The Hamiltonian 
of the system is 

i7 = J^S..S, = |^(Si-9/4), (3) 

where the Sa is the total spin of triangle A. The energy is minimized when Sa = 1/2 for every triangle, which can be achieved 
by putting a singlet bond on every triangle. The ground state is doubly degenerate. The two ground states shown in Fig. 2(b) and 
(c) violate the symmetry of reflection. 

Two types of domain walls interpolate between the ground states: the kink and the antikink. Fig. 2(d). A kink is an excitation 
with zero energy that happens to be an exact eigenstate of the Hamiltonian (3). Thus kinks are locahzed in the exchange-only 
model. The localized nature of kinks can be traced to an accidental degeneracy of the ground state of the exchange Hamiltonian 
on a triangle with half-integer spins in addition to the two-fold Kramers degeneracy. The two degenerate states with Sz=l/2 
have spin current going clockwise or counter clockwise around the triangle. The states also carry electric currents of opposite 
directions.'*^ An alternative set of basis states would have distinct valence-bond averages (S; • Sj) on the three bonds, which 
translates to nonzero electric charge on the three sites. '^^ 

In contrast, an antikink is mobile. The motion of an antikink is accompanied by the emission and absorption of kink-antikink 
pairs. The existence of a finite spin gap guarantees that these excitations are virtual. Polarization effects can be taken into 
account by using a variational approach. At the crudest level, the Hamiltonian (3) is projected onto the Hilbert space with a 
single antikink to obtain an effective hopping Hamiltonian for an antikink: 

K J J J 

H(^^\x) = -\x)~-\x + l)--\x-l). (4) 

where \x) is a state with an antikink on triangle x. The energy dispersion of the antikink is 

Ea{k) = 5J / 4 - J cos k, (5) 

with the minimum energy A = 0.25 J. In view of the zero energy of a kink, this value is the spin gap. 

This estimate can be further improved by enlarging the Hilbert space to include virtual excitations in the immediate neighbor- 
hood of an antikink. This yields an improved estimate of the spin gap, A = 0.219 J,"*" which is quite close to the result obtained 
by exact diagonalization, A — 0.215 J. "'^ 

It seems clear from the above that the variational approach provides a reliable description of the low-energy spin excitations 
in the pure Heisenberg model. We will use the lowest-order approximation for D ^ 0, without correcting for the vacuum 
polarization, to obtain a rough estimate for the critical coupling Dc- 

B. D/O 

In the presence of a nonzero DM term, kinks become mobile. For a single triangle, this means the splitting of the accidental 
degeneracy mentioned previously: the energy of a state with Sz = +1/2 now depends on the orbital momentum, reflecting the 
spin-orbit origin of the DM term. 

For an infinite chain, we follow the variational method described above and work in the Hilbert space spanned by states |.t) 
with a single kink located between triangles x and x + I. These states are not orthogonal to each other because they are not 
eigenstates of the same Hermitian operator The overlap is 

(xi|x2>= 2-1^1--^'. (6) 

As with antrkinks,'*" a simple rotation can be made to obtain an orthonormal basis {|i)}: 

|i) = A|:,)__l|^_l). (7) 

The matrix elements of the effective Hamiltonian in this subspace are 

{ii\H\i2) = -^2-1-1-^^1 sgn(xi - X2), (8) 



where the sign function is defined in such a way that sgn(O) = 0. A Fourier transform of the matrix element yields the energy 
dispersion of the kink: 

, , 6Dsiiifc 

The bottom of the band is at E"™'" = — 2|D|. For 13 > 0, it is reached for an incommensurate wavenumber k/2TT = 
-acos(4/5)/27r « -0.10. 

The calculation of the antikink case proceeds in a similar way. The basis states { |a;) }, with an antikink located at triangle x, 
can be orthogonahzed in the same way to yield an orthonormal basis { | £) } . The matrix element of the DM term is 



{S:i\Hdm\S:2) = -iD2-\^'-^^\ sgn{xi ~ X2) 
The resulting antikink dispersion is 



2 ~ o ('^3:i,X2 + l + (^xi.Ka-l) 



(10) 



r. n^ rr,. r , SZ? . , 313(4 COS fc - 1) sill fc 

E^{k)^5J'i-Jcosk + --smk+ '-— — i . (11) 

6 10 — 8 cos fc 

For D <^ J, the lowest energy of an antkink E™'™ = J/4 - lAD'^ / J + OiD"^ /J'^). The bottom of the band is located at 
fc/27r = -8D/3ttJ + 0{D^/J'^). 

The above energy dispersions were computed for spinons with Sz = +1/2. The dispersions for Sz = —1/2 can be obtained 
by changing fc i-> — fc. 

The bottom edge of the two-particle continuum as a function of total momentum is shown as solid lines in Fig. 3 for Sz ~ 
and in Fig. 4 for Sz = +1. (The former is a combination of two continua, one for a kink with 5^ = +1/2 and an antikink with 
Sz = —1/2, the other for a kink with Sz — —1/2 and an antikink with Sz — +1/2.) The edge dispersion mostly tracks the 
dispersion of the heavier particle, in this case the kink (9). The minimum energy of a kink-antikink pair 

^min ^ j^4 _ 2\D\ - 14L>V J + 0{D^/J^) (12) 

vanishes when the DM coupling reaches the critical strength Dc = 0.09,7. The total momentum of a Sz = +1 spinon pair 
with the lowest energy is fc/27r « —0.15. The gapless state arising at this critical point is expected to have transverse spin 
fluctuations with this wavenumber. The wavenumber of longitudinal spin fluctuations is determined by the bottom of the two- 
spinon continuum with Sz = 0, which occurs at fc/27r w ±0.06. 

III. EXACT DIAGONALIZATION 

To test the theory, we have performed an exact diagonalization study of the sawtooth chain with exchange and DM interactions. 
We worked with finite chains containing 2L sites in a system with L triangles with periodic boundary conditions. The length 
varied from L = 5 to 15. Both uniform and staggered DM interactions were investigated, with qualitatively similar results. 
Here we report on the uniform case only. For the largest system sizes, we employed the Lanczos algorithm, which provides 
convergent results for the ground state energy and a limited number of low-lying excitations. To reduce the size of the Hilbert 
space, we used the symmetry of translations along the chain and the 0(2) symmetry of spin rotations around the z-axis. 

Figure 3 shows the low-energy portions of the spectra in the Sz = sector for a chain with length L = 15 (30 sites), for several 
values of the DM coupling D. The invariance of the Hamiltonian (1) under time reversal symmetry (Sz >-^ —Sz, k -^ — fc) 
guarantees that the 5*^ = spectra are symmetric under mirror reflection (fc — > — fc). The lowest-energy excitations in the Sz — 
sector are expected to be spinon pairs in two channels: a kink with Sz = —1/2 and an antikink with Sz = +1/2 or vice versa. 
The calculated edges of the two-particle continua reproduce the shape of the dispersing bottom reasonably well. However, the 
calculated edge shifts downward with D faster than the numerical data do. 

In the Sz = +1 sector, the spectra are not symmetric under the mirror symmetry (the Sz = I spectrum maps onto that of the 
Sz = —I sector). Fig. 4. The lowest-energy excitations are expected to be spinon pairs consisting of a kink and an antikink, both 
with Sz = +1/2. Again, the calculated bottom edge of the excitation continuum has the right shape but advances downward 
with D somewhat too fast. In the two-spinon approximation, both the Sz ~ and Sz ~ I continua touch zero energy at 
Dc = 0.09 J. However, the numerical energy spectra appear to still have a gap at that point, see Fig. 3. 

To locate the critical point, we turned to a scaling analysis of the ground-state splitting. In the phase with valence-bond 
order, the ground state is doubly degenerate in the limit L ^^ oo. In finite systems, the ground-state doublet is split thanks to 
quantum tunneling. Both members of the doublet have momentum fc = because the valence-bond order preserves translational 
symmetry. The tunneling amplitude decays exponentially with the system length L and so does the splitting. 
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FIG. 3: Low-energy spectra of the sawtooth chain with a uniform DM term in the Sz = sector. Energy levels, measured relative to the 
ground state, are shown as a function of total momentum. Circles are the results of exact diagonalization for a periodic chain of length L = 15. 
Solid curves show the bottoms of the two-spinon continua computed analytically. Dashed straight lines show a linear dispersion with the speed 
V = 0.36J. 
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FIG. 4: Low-energy spectra in the Sz = +1 sector. Notations are the same as in Fig. 3. 



Fig. 5(a) shows the splitting of the ground state for D < O.llJ. All of the data sets, with the exception of the largest coupling, 
are well fit by the scaling expression 



AE = AL-^/^e-^/^ cos {kL) 



(13) 



with the same prefactor A. The dependence of the tunneling length ^ and the wavenumber k is shown in Fig. 5(c). The tunneling 
length diverges, or at least greatly exceeds the maximum attainable system length L = 15, for D > Dc = 0.115J. For 
O.llJ < D < 0.15 J, the finite-size dependence of the splitting was best fit by Eq. (13) with ^ = oo and a D-dependent 
amplitude A, Fig. 5(b). Apart from the oscillating factor, Eq. (13) suggests a scale-invariant ground state for D > Dc. The 
oscillations presumably come from the interference of instantons as discussed in the Appendix. 
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FIG. 5 : The splitting of the ground-state doublet as a function of the system length L for (a) D < Dc = 0. 1 15 J and (b) for D > Dc ■ (c) The 
dependence of the inverse tunneling length 1/^ and the wavenumber k in the scaling form (13) on the DM coupling strength D. 

For D > Dc, we expect a gapless phase with quasi-long-range incommensurate spin correlations decaying as a power of the 
distance. For a sufficiently large D, the classical model should become a good starting point. In the classical limit, the sawtooth 
chain has a spiral order for any nonzero value of D, Fig. 2(f). Low-energy excitations are spin waves with a speed 



2.7S\ 



(14) 



Quantum fluctuations disrupt the long-range spin order, restoring translational invariance and the 0(2) symmetry. Such a phase 
would be a Luttinger liquid, whose lowest-energy 5*2 = +1 excitations are spin waves with a sound-like spectrum at fco/27r = 
— 1/3. The numerically determined Sz = +1 spectra for D > 0.15J are consistent with spin waves. At D = 0.19J, the soft 
spot is located at fco/27r « —0.25, not far from the classical value. The speed of sound (estimated from the slope of the dashed 
lines in Fig. 3 and 4) is w = 0.36 J, is not far from the classical estimate (14) obtained below. 

IV. SPIN CORRELATIONS IN THE GROUND STATE 

To verify the location of the quantum critical point Dc and to confirm the critical nature of the ground state for D > Dc, we 
examined the long-distance behavior of spin correlations, G'^^{r) ~ {S"{0)S^{r)), in the ground state. In the Luttinger-liquid 
regime, transverse spin correlations are expected to decay as a power of the distance,"*"^ 

^MI-OT- (15) 



\G^ 



The stiffness constant K varies between 1 (gas of dilute magnons) and 1/4 (gas of dilute spinous). ^^'^^ 

In a finite system of length L with periodic boundary conditions, the Green's function depends in the same way on the chord 
distance"*^ 



d(r) = (L/tt) sin (7rr/L). 



(16) 



In a system with 2L spins, this distance varies from rf w 1 to L/tt. In view of that, the range of distances in a system with 
2L = 30 spins is not sufficient to reliably observe the critical behavior of the spin correlation function. 

To observe the critical behavior, we used the density-matrix renormalization group (DMRG) method implemented through the 
Matrix Product Toolkit"*^ to obtain the ground-state wavefunction in a periodic chain with up to 2L = 100 spins. The system has 
a U(l) symmetry which we took into account to reduce CPU time. The number m of states kept varied from 800 to 1200 states. 
Our results for the ground state energy per site for all values of DM coupling D investigated are consistent with the energy per 
site obtained from the ED calculations. 

The resulting transverse spin correlations jG^ (r)| in a system of length L = 50 are shown in Fig. 6 as a function of the 
chord distance (16). At largest distances d, the data for D = 0.12 J follow a power law C /<P, which is consistent with the 
value K = 1/4 at the spinon condensation point. For D > 0.12J, spin correlations follow power laws with smaller slopes, 
indicating K > 1/4. For D < 0.12 J, the power-law scaling breaks down at large d changing to an exponential dependence. 
The estimated critical point, Dc = 0.12 J, is in reasonable agreement with the value Dc = 0.1 15J obtained from the splitting of 
the ground-state doublet. 



V. DISCUSSION 



Analytical arguments and numerical evidence presented above supports the following scenario. In the absence of the 
Dzyaloshinskii-Moriya term, the sawtooth chain has a doubly degenerate ground state with valence-bond order spontaneously 
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FIG. 6: The amplitude of transverse spin correlations ( 1 5) as a function of the chord distance ( 1 6) on a log-log plot (left) and a simple log plot 
(right). 



breaking the reflection symmetry of the lattice. Elementary excitations are spinons of two flavors, localized kinks and mobile 
antikinks. The gap to spin-1 excitations, A = 0.215 J is determined by the edge of the two-spinon continuum. The introduction 
of a DM term with the D vector pointing along the same axis for all bonds. Fig. 2, lowers the spin-rotation symmetry down to an 
0(2). At weak coupling D, the lattice reflection symmetry remains spontaneously broken. At the same time, a finite D lowers 
the excitation energies of both kinks and antikinks and the spin gap (understood as the lowest energy of 5"^ = 1 excitations) 
begins to close. A fairly crude analytical calculation indicates that the main factor affecting the spin gap is the minimum energy 
of the kink, — 2|D|. The gap closes roughly when that energy equals the initial gap in absolute terms, \D\ = Dc ^ A/2 w O.IJ. 
This is confirmed by numerical work involving exact diagonalization of finite chains, with the result Dc — 0.115 J. Beyond the 
critical coupling, the spinons proliferate. Since they act as domain walls in the valence-bond order parameter, the valence-bond 
order is lost and the lattice symmetry is fully restored. The resulting state is likely a Luttinger liquid with incommensurate spin 
correlations and spin-wave excitations. Similar transitions between Ising-ordered phases and Luttinger liquids have been found 
in other one-dimensional systems.'*''*^ The strength of the DM coupling D w {5g/g) J where Sg is the deviation of gyromagnetic 
ratio from its free-electron value g.'^ In kagome antiferromagnets herbertsmithite and volborthite, 5g/g w 0.1."*^ 

It is tempting to speculate that a somewhat similar transition may occur in the 5* = 1/2 Heisenberg model on kagome with 
a DM coupling. While the existence of the transition is not in doubt — at a large enough D the system should develop magnetic 
order^^"^^ — the nature of the transition remains to be determined. 

In the kagome antiferromagnet, spinon excitations are very similar to those of the sawtooth chain."*" In the absence of the DM 
term, kinks are localized and have zero energy, whereas antikinks follow one-dimensional trajectories with the same energetics 
as on the sawtooth chain. Adding the DM term thus has similar consequences, namely delocalization of kinks is the main factor 
lowering the edge of the kink-antikink continuum. If anything, the gap may close even faster than on the sawtooth chain because 
on kagome kinks move in two dimensions and thus can lower their energy through delocalization more effectively than on a 
one-dimensional chain. For this reason, the critical DM coupling for kagome may be even lower than for the sawtooth chain. 

The kagome antiferromagnet differs from the sawtooth chain in one important respect: it has a finite concentration of antikinks 
in the ground state. The antikinks form tightly bound S = pairs, whose binding energy Aaa ~ 0.06J is lower than the 
threshold energy of kink-antikink creation Aka ~ 0.25 J. Therefore the spin gap in the Heisenberg antiferromagnet on kagome 
is determined by binding energy of an antikink pair Although the binding energy Aaa is no doubt influenced by the introduction 
of the DM term, it is unlikely that this energy is very sensitive to the presence of a small perturbtion like D as Aaa is determined 
by a competition of two high-energy processes: the antikink hopping amplitude and the antikink attraction in the singlet channel, 
both with a strength of order J. It seems more likely that the larger gap Aka will be quickly driven to zero as it is on the sawtooth 
c hain. 

The nature of the phase transition at the conjectured condensation of kinks and antikinks is an open question. It is not even 
known whether the D = ground state is a valence-bond liquid or solid, with contradictory indications from different numerical 
techniques.^'"-"* (In our view, even a small amount of bond disorder will turn the system into a disordered valence-bond solid.) 
Adding the DM term will tend to melt the delicate valence-bond order turning the valence-bond crystal into a liquid before 
the magnetic condensation and thus inducing another phase transition along the way. The nature of the condensed phase is 
not clear, either Usually, ordering of the transverse components of magnetization is associated with a proliferation of 5"^ = 1 
objects, as is the case in magnon condensation,-"' whereas here the condensing particles are spinons with half-integer spin. This 
obsrvation lends support to the scenario with an intermediate gapless phase lacking long-range spin order, ■'^ which is some sort 
of an algebraic spin liquid."* 
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Appendix A: Oscillations in the ground-state splitting 

To understand the oscillatory behavior of the ground-state splitting, Eq. (13), we turn to a much simpler model: the antiferro- 
magnetic XXZ chain with DM interaction described by the Hamiltonian H = -ffxxz + ^^dm, where 

Hxxz = E [JcosaiS:s:+, + SyX^,) + J.S^^S-^^,] (Al) 

n 

and 

iJoM = JsmaJ2iS^S'n+, - ^,".^,"+i). (A2) 

n 

In the easy-axis limit, J^ S> J, the ground state is doubly degenerate and exhibits Neel order. In a finite chain with periodic 
boundary condition, quantum tunneling splits the doublet into eigenstates with momenta and tt. Below we discuss the effect 
of the DM term, a 7^ 0, on the splitting. 

By rotating local axes at site n through angle na in the xy plane, the DM term in the Hamiltonian can be removed, producing 
the standard XXZ model: 

H = 2_^[J\^n^n+l ^ ^n^n+l) ^ 'Jziin^n+li ■ (^3) 

n 

For a closed chain of length L, the transformation yields twisted periodic boundary conditions: 

5"^ = ^o cos (La) + 5"^ sin (La), 

S% = - S^ sin (La) + S^ cos (La). (A4) 

The twist is absent if La = 27rm, where m is an integer. Then the system has the same spectrum as in the absence of the DM 
term, a = 0. At a fixed chain length L, the splitting is a periodic function of a with a period of 2tt/L. 

To see that the splitting should have an oscillatory character, consider the special case of a tt twist. La = {2m + 1)tt. As 
Haldane argued,^^ the tunneling between the two Neel states is mediated by instantons with quantized winding numbers n, 
classical action S'„, and a Berry phase cxp (27rm5). For boundary conditions with a tt twist, the winding numbers are half- 
integer, 77, = ±1/2, ±3/2, . . . Instantons with opposite winding numbers have the same classical action, Sn = 5_„. However, 
their Berry phases are exactly opposite, cxp {2'KinS) = — cxp (— 27rm5), when both the winding numbers n and spin S are half- 
integer As a result of destructive interference of instantons with opposite winding numbers, the tunneling amplitude vanishes 
when La = {2m + l)7r. We thus expect an oscillatory dependence of the splitting on a at a constant L in the XXZ chain with 
half-integer spins and periodic b oundary conditions. The exponential dependence of the splitting on the length will acquire an 
oscillating prefactor cos {aL). This inspired Eq. (13). 

Appendix B: Spin wave in sawtooth chain 

We compute the spin-wave spectrum on the sawtooth chain in the classical limit, 5 ^ cx). The Hamiltonian is 

i/ = ^[S,.S,+D„-.(S, xS,)]. (Bl) 

For brevity, we set J = 1. 

In equilibrium, spins lie in the plane normal to the DM vectors D^, with the angle of 120° between nearest neighbors. 
Fig. 2(f). It is convenient to choose reference frames in such a way that spins point along the local z axes, the x axes are in the 
plane of the spins, and the y axes are parallel to D^ . For small deviations from equilibrium, 

S, « S{a,, /3„ 1 - ajl2 - p^j2) (B2) 



10 

where a^ and /?i are small deviations from the 120° pattern. 
In the harmonic approximation, the energy (B 1 ) reads 

where x = ~l/2 — V3D/2. Ki = 1 if z is an apex (A) site and Ki = 2\i it is a base (B) site. 

The dynamics can be obtained from the Lagrangian, which includes a Berry phase term in addition to the potential energy: 

L = sY{cose,-l)^,-H. (B4) 

i 

After expressing the angles 9 and in terms of a and /3, 

tan(/) = /?/a, cose* « 1 - (a^ + /3^)/2, (B5) 

we obtain the following Lagrangian: 

L = sY,{aA ~ a^A)/2 - H. (B6) 

i 

It yields the equations of motion for spins on sublattices A and B\ 

at = 5x(/3.>i/2 + /^f-i/2)-2%/3f, (B7a) 

af = SxiPti/2 + Pti/2 + PF+1 + /3f-i) - 45x/3f , (B7b) 

Pt = -5«i/2 + af_i/2) + 25x«f , (B7c) 

/3f = -5(af+i/2 + af_i/2 + <i + «f_i) + 45xaf . (B7d) 

Note that i is half-integer on sublattice A and integer on sublattice B. Plane waves with frequency a; and wavevector k satisfy 
the equation 



~iuj 2Sx -2Sxcos{k/2) \ f a^ \ 

-iuj ~2Sxcos{k/2) 4S'x - 2S'x cos fc 

-2S'x 25008(^/2) -iuj 

2Scos{k/2) -ASx + 2S cos k -iuj 






(B8) 



At _D = 0, we have one zero mode and one mode with a finite frequency w = 5*^2 — cos(2fc). For a finite D, the zero mode 
acquires a dispersion linear in k in the limit fc ^- 0. The wave velocity is 



V 2 + 8\/3L> + 18L»2 

Restoring J as a coupling constant, we find the following behavior for the velocity. As Z? ^ 0, u ^ 2.79S^/DJ. For 
D = 0.19J,v = 1.05SJ. 
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